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The extraction of scattering parameters from Euclidean simulations of a Yukawa model 
in a finite volume with periodic boundary conditions is analyzed both in non relativistic 
! quantum mechanics and in quantum field theory. 
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1. Introduction 

A rigorous treatment of nuclei based on first principles quantum chromodynamics re- 
■ quires the understanding of the transition from the confining theory of quarks to the 
t> . unconfined hadrons, a present challenge for the high energy community. Only very re- 
cently, the development of algorithms and computers devoted to lattice QCD has made 
possible starting the study of nucleon-nucleon [1] systems (as well as hit, KK, tcN, etc) 
based on full QCD calculations. Nevertheless, lattice methods have a drawback for com- 
puting scattering parameters, as they can not be extracted in the infinite volume limit [ 
|2j. However they can be obtained by measuring the finite volume effects on the spectrum 
of a periodic boundary condition box 

In this work, the scattering length of the Yukawa model has been determined; on 
one hand, for the Yukawa potential V(r) = —ae~ msr /r in Schrodinger equation, that 
can be precisely solved in a 3-d torus, as well as in the infinite volume, serving as a 
crosscheck for determining scattering length via finite volume effects. On the other hand, 
the quantum field theory of fermions (ip) and mesons ((f)) interacting via the lagrangian 
Cj(x) = gQip(x)Y(j)(x)il)(x) - that gives rise to the latter potential (with a = gljAm) in the 
NR limit - is solved in a Euclidean lattice. 

The goals of this calculation are manyfold. First, the comparison between the QFT and 
the non relativistic (NR) results can help to disentangle the bias introduced by the ladder 
approximation, as well as the size of the relativistic corrections. Second, the methods used 
to extract the scattering length can be tested with this model, much simpler than QCD. 
And finally, the interaction itself is interesting, as scalar meson exchange is a common 
ingredient of NN interaction potentials [ 0] , responsible for most of the nuclear binding 
energy. 
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2. Prom bound state spectrum to low energy parameters 



The relation between the spectrum of a system enclosed in a torus and the phase-shifts 
in the infinite volume has been established by M. Luscher [EJ]. For a lattice with L points 
in each spacial direction and spacing a, the energies are determined by the expression: 
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that is exact for lattice sizes aL > 2R where R is the interaction range, i.e. V(x > R) = 0. 
For interactions of physical interest this regime is reached exponentially and independently 
of the coupling constant and the low energy parameters (LEP) appearing in the effective 
range expansion, pcot(8(p)) ~ — + \ r oP 2 + ■ • •■ This constitutes a remarkable advantage 
with respect to the, most commongly used, large-L expansion of (jTj): 
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G=0.4, am =0.2 



In a recent work the possibility to ex- 
tract the LEP from the bound state spec- 
tra in a finite box has been examined in the 
framework of the NR Yukawa model. The 
L-dependence of the eigenenergies E(L) in a 
3d torus was used to extract the infinite vol- 
ume scattering length ao and effective range 
r . Contrary to the quantum field lattice 
calculation, those quantities can be indepen- 
dently computed by solving the correspond- 
ing Schrodinger radial equation and used to 
test the validity of the different approaches 
discussed in the literature. 

The results were analyzed in terms of 
the slowly varying quantity Oq\l) = 
(aL) 3 E(L)M/4:7r which represents the zeroth- 
order approximation of large-L expansion (J2J). 
The NR model depends on the unique param- 
eter G = aM/m 
- according to (J2J) 
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Figure 1: Lattice result for <Zq vs volume 
compared to the infinite volume result. 



s , and it was shown that a 
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(L)m s tends to —G in the small-L limit, and 
to the required scattering length value a for aL — > oo (figure [TJ). It 
is interesting to remark that a^ (L) displays the same 1/L 3 behavior in both limits. 

It was in particular found that one can obtain accurate values of ao and ro by using 
equation (JIJ at two different lattice sizes in the region Lam s w 5 and solving the resulting 
linear system. Some results are illustrated in figure El They show the sensibility to the 
effective range ro and the need for both parameters to be simultaneously determined. 

The possibility of extracting LEP that way is independent of the ao value and applies 
in the resonant well. Care must be taken however when applying equation (J2J to 

extract the value of a . When using lattice sizes Lam s ~ 2—3, E(L) is well fitted by a c/ L 3 
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dependence but the coefficient c can strongly differ from the infinite volume scattering 
length Only the use of equation (JTJ) at lattice sizes greater than the interaction 

range - Lam s > 5 in the Yukawa model - could lead to unambiguous extraction of LEP, 
provided both a and r are taken into account. 



3. Lattice results 

The non perturbative solutions of the 
Yukawa field theory in the lattice have been 
considered in [ El El • The existence of a criti- 
cal value of the coupling has been argued [ 
IB] hindering the appearance of bound 
states of the QFT. The effect of the ladder 
approximation can nevertheless be studied by 
calculating the scattering length and compar- 
ing with the non relativistic results. 

We have used Wilson discretization of 
fermion fields, which has the disadvantage 
that there are 0(a) discretization errors. 
They can be important in our case, since nu- 
cleon mass is not small. As the interaction 
is non confining, to reduce the discretization 
errors we will use a hint from the free case, 
where the time evolution of fermion is given 
by 



-10 - 



-15 ■ 



-20 



-25„ L 



G=1.5, am =0.10 



•- -. r fitted 

a- a r„m =2.46 (exact) 

— Infinite volume 



10 

aL 



15 



20 



Figure 2: ao obtained with equation £7p for 
different values of ro . The influence of this 
parameter is remarkable. 



C(t) = ^2Tr(ip(x,t)'ip{d,0)) oc(l + aM L )"'~ e 
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This expression defines the mass that governs the time evolution of the nucleon M, and 
relates it to the fermion mass appearing in the lagrangian, Ml, as aM^ = exp(aM) — 1. 

Both masses are equivalent at the continuum limit, but they can differ sizably when 
working at finite lattice spacing, what may spoil the convergence of the lattice results. 
To overcome this, we define an improved coupling, g, as the variation of nucleon energy 
with the meson field at zero meson field (note that the coupling acts as a mass term), 
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If we assume that the nucleon propagates in a small constant background field, fy, its 
mass is modified by the coupling term as aM^ —>■ aM^ + g o 0& which implies aM — > 
aM + g (f>b/(l + aMj)j. The renormalized coupling defined above results then 
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Note that go an d g are equivalent in the continuum limit (a — > 0), but the second should 
have a better behavior at finite lattice spacing. The usefulness of this coupling is specially 
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important when working with coarse lattices, as is the case in nuclear physics, due to the 
mass hierarchy. Note that the differences due to this definition of the coupling can be 
very important, ~ 70% for fermion masses aM ~ 0.3. 

Two body S-wave scattering length can be 
computed from the corresponding spectrum 
on a finite size lattice calculation as it has 
been shown in the NR potential model. The 
fist test to be done is the perturbative regime 
that, as in the NR case, predicts for the scat- 
tering length aofrig = —G + 0(G 2 ). The re- ^S" 
suits of the scattering length obtained from "~«s 
lattice calculations are shown in figure El as a 
function of the coupling constant a = g 2 /4:ii, 
both for the original and the improved cou- 
pling. The correspondence with the small 
coupling limit (dashed line) gets clearly bet- 
ter with the use of the improvement factor 
®. 

The fact that we do not see differences Figure 3: Scattering length versus G for 
with the first order in perturbation theory the original g and modified coupling g for 
is maybe due to the fact that we work with lattice volumes Lam s ~ 2 - 4. 
rather reduced lattice volumes. Further cal- 
culations have to be done (in particular bigger volumes with La > 5m £ T 1 ) in order to 
compute the scattering length - and effective range - properly. 

As a concluding remark, we have proposed an improved definition of Yukawa coupling 
constant in the lattice that exhibits quite better scaling properties than the bare one. 
As shown in figure El this correction factor is mandatory to test perturbation theory - 
and accordingly to obtain any physical result beyond this regime - even for rather small 
lattice spacings aM ~ 0.2. 
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